Statistical and dynamical fluctuations in the ratios of higher net-proton cumulants 

in relativistic heavy ion collisions 
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With the help of transport and statistical models, we find that the ratios of higher net-proton 
cumulants measured at RHIC are dominated by the statistical fluctuations. Future measurements 
should focus on the dynamical fluctuations, which are relevant to the underlying mechanisms of 
particle production, the critical phenomena in particular. We also demonstrate that a proton- 
antiproton correlation directly show if protons and antiprotons are emitted independently. 
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I. INTRODUCTION 

One of the main goals of relativistic heavy ion collisions 
is to locate the critical point on the QCD phase diagram, 
spanned by the temperature (T) and the baryon chemi- 
cal potential (/is). At the critical point, the correlation 
length (£) goes to infinity and the long range correlations 
become dominant. So the fluctuations of final state par- 
ticles are expected to be largely enhanced, if the (T, he) 
trajectory of the collision system is close to the critical 
point. The ^-related observables are therefore of great 
interest in heavy ion collisions 

The thermodynamic quantities, such as order parame- 
ter, specific heat capacity and susceptibility (x), diverge 
with the correlation length at the critical point. The 
ith net-baryon cumulant is recently shown to be directly 
related to the ith susceptibility of the formed sys- 
tem dm, 



(5N l ) = VT Xi , 



(1) 



where N is the net-baryon number, (6N l ) = ((N—(N)) i ) 
is the ith net-baryon cumulant, and V is the volume. The 
third and fourth cumulants, 



A 3 = (SN 3 ), K 4 = (5N 4 ) - 3(5 N 2 ) 2 , 



(2) 



arc argued to be more sensitive to the correlation length 
as they are proportional to £ 4 ' 5 and £ 7 , respectively p-Q- 
Experimentally, the proton number is a good approxima- 
tion of the baryon number Q, and the properly normal- 
ized ratios, net-proton Skewness and Kurtosis, 
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are preferred Q , measuring the symmetry and sharpness 
of the net-proton distribution, respectively. 

The STAR measurements Q show that both net- 
proton Skewness and Kurtosis decrease with increasing 



number of participants (centrality), which could be ex- 
plained by the central limit theoremfCLT). On the other 
hand, various model calculations Q reproduce the ex- 
perimental results surprisingly well, raising the suspicion 
that Skewness and Kurtosis are insensitive to the mech- 
anisms of particle production implemented in different 
models. Recently, Karsch and Redlich Q have derived 
simple relations between the cumulant ratios and the 
thermal parameters (T and /is at the chemical freeze- 
out flCI). based on the hadron resonance gas (HRG) 
model [ill ], with the system well thermalizcd and with- 
out phase transition. They have shown that the HRG 
model results are well consistent with the STAR data at 
different collision energies @. 

However, before trying to understand the physics be- 
hind the measured Skewness and Kurtosis, we need to 
take into account and properly remove the contribu- 
tions from the statistical fluctuations [!, [l2j and all non- 
thermal sources (miniiets, resonance decay, initial size 
fluctuation, etc.) [HI, Q|- I n this paper, we will focus 
on the elimination of the statistical fluctuations. We 
first estimate the statistical fluctuations in the cumu- 
lant ratios, which turn out to dominate the behavior of 
net-proton Kurtosis at RHIC energies. Then we propose 
the dynamical ratios of higher net-proton cumulants in 
Section III and the correlation between proton and an- 
tiproton in Section IV. The centrality dependence of the 
dynamical ratios and the correlations from two versions 
of AMPT [3, UrQMD 0, and Therminator [l7| are 
presented and discussed. Finally, the summary and con- 
clusions arc given in Section V. 



II. STATISTICAL FLUCTUATIONS IN THE 
RATIOS OF HIGHER NET-PROTON 
CUMULANTS 



The statistical fluctuation comes from the finite num- 
ber of particles, usually obeying a Poisson distribu- 
tion 0, [L2|, [l9| . If we have two independent Poisson dis- 
tributions for protons and antiprotons with means (N p ) 



2 



and (Np), respectively, the net-proton number (N) fol- 
lows a Skellam (SK) distribution HGJ], 

f(N; (N p ),(N p )) 

= e ~((N p)+{ N P )) {{Np)/{Np)) f IiNl (2yJ(N p )(N p f) , 

(4) 



where Lj^i (2-^/ (N p ) (Np) ) is the modified Bcssel function 
of the first kind. Then the statistical fluctuations of the 
ratios of higher net-proton cumulants can be directly de- 
duced from the Skellam distribution, 
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where Rij = Xi/Xj — Ki/Kj- The HRG model re- 
sults Q show that Z^4 t 2 is unity, and now we find that it 
is completely statistical. The other ratios are determined 
only by the mean numbers of protons and antiprotons, 
which usually increase with the incident energy and cen- 
trality. So the statistical fluctuations of Skcwness and 
Kurtosis decrease with the incident energy and central- 
ity. 

In Fig. 1(a), the model results of net-proton Kurto- 
sis and its statistical fluctuations (SK) are shown from 
APMT default, AMPT with string melting [H[ and Ther- 
minator (l7| . Where the transverse momentum, p t , and 
rapidity cuts are respectively 0.4 < p t < 0.8 GeV/c, 
\y\ < 0.5, the same as they are given at RHIC/STAR 
published paper 0] . In both transport (AMPT) and sta- 
tistical (Thcrminator) models, the net-proton Kurtosis 
results (solid points) are very close to the corresponding 
statistical fluctuations (open points). The STAR mea- 
surements of net-proton Kurtosis (stars) Q follow the 
same trend as a function of the number of participants 
(N part )- So the statistical fluctuations dominate the be- 
havior of net-proton Kurtosis at RHIC, which is why the 
results from various models closely resemble the experi- 
mental data. 

In Fig. 1(b), we show the incident-energy dependence 
of the net-proton ratios, R.2,1, Ri,2, and i?3,2, obtained 
from RHIC/STAR data, HRG, AMPT default, AMPT 
with string melting and Thcrminator, together with the 
corresponding statistical fluctuations (SK) of the last 
three cases determined by Eq. (5). At 200 GeV, the ra- 
tios from Thcrminator (squares) coincide with the corre- 
sponding statistical fluctuations (crosses), because Ther- 
minator is a statistical model with only constraints on 
kinetics. The HRG curves @ are given by i?2,i = 
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FIG. 1: (Color online) (a) Centrality dependence of net- 
proton Kurtosis at 200 GeV, and (b) Energy dependence of 
R2,i, Ra,2 and Rs,2, for Au + Au collisions. The results are 
obtained from RHIC/STAR data 0], HRG [jj, two versions 
of AMPT [l5|j and Therminator [l7| models, and the corre- 
sponding statistical fluctuations are estimated with the Skel- 
lam distribution (SK) from Eq. {^J, respectively. 



1/ tanh(/xs/T) and R3.2 = tanh(/is/T). Similarly, both 
Thcrminator and HRG models have protons and antipro- 
tons completely independently emitted with Poisson dis- 
tributions, so their results are close to the data from 
RHIC/STAR, where the statistical fluctuations domi- 
nate. The ratios from two versions of AMPT (solid cir- 
cles and triangles) slightly deviate from the correspond- 
ing statistical fluctuations (open circles and rhombi) , and 
from the HRG lines, due to the more complicated particle 
production mechanisms in AMPT. 

We have demonstrated that the influence of statistical 
fluctuations is far from being negligible in the ratios of 
higher net-proton cumulants at RHIC energies. To in- 
vestigate the underlying physics, we have to first remove 
the statistical fluctuations. 



III. DYNAMICAL RATIOS OF HIGHER 
NET-PROTON CUMULANTS 

There have been long efforts in eliminating the statis- 
tical fluctuations in elementary collisions [TH, For 
a single particle distribution, the factorial moments are 
used to remove the statistical fluctuations 0, [ljj • But 
this method can not be directly generalized to the distri- 
bution of the difference between two Poisson variables. 

From previous discussions, the statistical fluctuations 
in the ratios of net-proton cumulants are directly obtain- 
able. The dynamical ratios of net-proton cumulants can 
be simply defined as a deviation of the ratios from the 
statistical fluctuations [l9|], e.g., 



K, 



dyn 
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(6) 



and so on, where the statistical parts are given by Eq. . 

The centrality dependence of dynamical net-proton 
Skewness and Kurtosis are shown in Fig. 2(a) and (b), re- 
spectively, from AMPT default, AMPT with string melt- 
ing, UrQMD and Therminator models for Au + Au colli- 
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FIG. 2: (Color online) Centrality dependence of dynamical 
Skewness (left) and Kurtosis (right) for Au + Au collisions at 
200 GeV, given by transport models (AMPT and UrQMD) 
and a statistical model (Therminator) . 
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FIG. 3: (Color online) Centrality dependence of dynamical 
net-proton Skewness (left) and Kurtosis (right) for Au + Au 
collisions at 3 incident energies, obtained from AMPT default. 



sions at 200 GeV. Both dynamical Skewness and Kurtosis 
from Therminator are zero at all centralities, illustrating 
that the symmetry and sharpness of the net-proton distri- 
bution in the model both follow the Skcllam distribution. 
For the transport models (AMPT and UrQMD), both 
dynamical Skewness and Kurtosis are larger than zero in 
peripheral collisions, and approach zero in central colli- 
sions. Compared with the Skcllam distribution, the pos- 
itive dynamical Kurtosis and Skewness implies that the 
net-proton distribution has a sharper peak and a longer 
tail at the large net-proton side, respectively. These de- 
viations are caused by non-thermal sources implemented 
in transport models. 

To study how the results change with the incident en- 
ergy, the centrality dependence of dynamical net-proton 
Skewness and Kurtosis from AMPT default are shown in 
Fig. 3 for Au+Au collisions at 3 incident energies. When 
the incident energy changes from 200 GeV to 39 GeV, 
both dynamical Skewness and Kurtosis remain positive. 
Dynamical Skewness shows a significant dependence on 
the incident energy, especially in peripheral collisions, 
and dynamical Kurtosis is almost independent of the in- 
cident energy. 

The behavior of dynamical ratios of higher cumulants 
are called for at the RHIC beam energy scan. If the 
deviation from the statistical fluctuations is zero, pro- 
tons and antiprotons are emitted independently as the 
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FIG. 4: (Color online) Centrality dependence of proton an- 
tiproton correlations for Au + Au collisions at 200 GeV, given 
by transport models (AMPT and UrQMD) and a statistical 
model (Therminator). 



statistical models assume. Otherwise if the non-zero de- 
viation remains the same sign for different incident en- 
ergies, like what the transport model shows in Fig. 3, 
then there is no critical related phenomena. However, 
if the deviation changes dramatically with the variation 
of the incident energy, e.g. showing the sign changes at 
the third and fourth cumulants, where the symmetry and 
sharpness of the net-proton distribution deviate from the 
corresponding statistical fluctuations in opposite direc- 
tions I6l. l2ll . |22j , it may reveal the critical incident energy 
nearby |22j. 



IV. CORRELATIONS BETWEEN PROTON 
AND ANTIPROTON 



To see if protons and antiprotons are emitted inde- 
pendently, we could also directly measure the correlation 
between them, 



C(N P ,N P ) 



(N P )(N p ) 



(7) 



The correlation will be zero if protons and antiprotons 
are independent. 

The centrality dependence of the correlations from two 
versions of AMPT, UrQMD and Therminator models for 
Au + Au collisions at 200 GeV are presented in Fig. 4. 
The correlation is zero at all centralities in Thermina- 
tor, another illustration of the model's assumption. In 
transport models, the correlation decreases with central- 
ity, following a similar trend as dynamical Skewness or 
Kurtosis. The correlation is positive, indicating that pro- 
tons and antiprotons are not emitted independently. This 
leads to the difference between the net-proton distribu- 
tion in transport models and the pure statistical Skellam 
distribution. 
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V. SUMMARY 

We have demonstrated that the statistical fluctuations 
dominate the behavior of the ratios of higher net-proton 
cumulants measured at RHIC, and this explains why the 
results from various models are consistent with the exper- 
imental data. We argue that before trying to understand 
the underlying physics the statistical fluctuations should 
be taken into account. 

To study the particle production mechanism, the dy- 
namical ratios of higher net-proton cumulants and the 
correlations between proton and antiproton have been 
proposed and discussed. It is shown that the dynamical 
ratios and the correlations are similarly zero at all cen- 
tralitics in a statistical model, and positive in transport 
models. This indicates that protons and antiprotons are 



not emitted independently in transport models. 

The behaviors of the dynamical ratios of higher 
net-proton cumulants, as well as that of the proton- 
antiproton correlation, are more relevant to the location 
of the critical point, and the corresponding measurements 
during RHIC beam energy scan will shed light on the 
study of the QCD phase transition. 
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